On the diffusive anomalies in a long-range Hamiltonian system 
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We scrutinize the anomalies in diffusion observed in an extended long-range system of classical 
rotors, the HMF model. Under suitable preparation, the system falls into long-lived quasi-stationary 
states presenting super-diffusion of rotor phases. We investigate the diffusive motion of phases by 
monitoring the evolution of their probability density function for large system sizes. These densities 
are shown to be of the g-Gaussian form, P(x) oc (1 + (q — l)[x/ P] 2 ) 1 ^ 1 ^^, with parameter q 
increasing with time before reaching a steady value q ~ 3/2. From this perspective, we also discuss 
the relaxation to equilibrium and show that diffusive motion in quasi-stationary trajectories strongly 
depends on system size. 
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I. INTRODUCTION 

Systems with long-range interactions constitute a very 
appealing subject of research as they display a variety 
of dynamic and thermodynamic features very different 
from those of short-range systems treated in the text- 
books (see [if for a review on the subject). Moreover, in 
recent years, the study of long-range models have raised 
a renewed interest due to the possible applicability of 
"Nonextensive Statistics" 2] to such systems. 

A very simple model that offers the possibility of inves- 
tigating many issues related to long-range interactions is 
the Hamiltonian Mean-Field (HMF) model 3]. It consists 
of N planar classical spins interacting through infinite- 
range couplings. The dynamical variables of each spin 
i are a phase angle 9i and its conjugate momentum pi 
whose evolution derives from the Hamiltonian 

N N 
i=l i,i=l 

This model can be seen as a variant of the XY ferro- 
magnet, equipped with a natural Newtonian dynamics. 
Although the range of interactions is infinite, the HMF 
has been shown to behave in many aspects qualitatively 
like its long(finite)-range analogs 4]. Then, despite its 
simplicity, it reflects features of real systems with long- 
range forces such as galaxies and charged plasmas 0. 

Its equilibrium thermodynamics can be solved in the 
canonical ensemble. It presents a ferromagnetic second- 
order phase transition, from a low energy clustered phase 
to a high energy homogeneous one that occurs at criti- 
cal temperature T c = 0.5 (i.e., critical specific energy 
e c = 0.75)0. However, concerning the microscopic dy- 
namics, the system may get trapped into trajectories over 
which averaged quantities remain constant during long 
periods of time with values different from those expected 
at equilibrium. For instance, for water-bag initial condi- 
tions (i.e., 0i — Vz, and pi randomly taken from a uni- 
form distribution), a quasi-stationary (QS) state appears 



at energies close below e c @. In a QS state, the temper- 
ature (twice the specific mean kinetic energy) is almost 
constant in time and lower than the canonical value to 
which it eventually relaxes. However, the duration of QS 
states increases with the system size N, indicating that 
these states are indeed relevant in the (N — > oo) thermo- 
dynamical limit (TL). 

Several other peculiar features have been found for 
out-of-equilibrium initial conditions, e.g., negative spe- 
cific heat[3j, non-Maxwellian momentum distributions |5j, 
glassy dynamics @, aging 0, anomalous diffusion 0, 
and others. In particular, anomalous diffusion has been 
initially associated to quasi-stationarityj9j and later to 
the (non-stationary) relaxation to equilibrium More- 
over, controversies about the characterization of these 
states have arisen[n|, remaining still aspects of the 
anomaly in diffusion to be investigated. In this work we 
report new results on anomalous diffusion and relaxation 
to equilibrium, focusing on the dependence on system 
size. 

II. RESULTS 

The equations of motion derived from Eq. 

Q\ = pi , for 1 < i < N 

Pi = My cos l - M x sin i} (2) 

where M = jf^2j{cos8j,sm9j) is the magnetization, 
were solved by means of a symplectic fourth order al- 
gorithm [T^ . Integrations were performed for fixed 
e = 0.69, for which quasi-stationary effects are more pro- 
nounced. Moreover, in the continuum limit, QS trajec- 
tories at energies around that value are stable station- 
ary solutions of Vlasov equation^! Q|. We considered 
two classes of initial conditions. One of them is a slight 
variation of water-bag initial conditions normally used in 
the literature: setting 6i = for i = 1, N, and regu- 
larly valued momenta (instead of random ones) with the 
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addition of a small noise to allow for statistical realiza- 
tions. This type of initial condition preserves the main 
features of random water-bag ones but the QS tempera- 
ture Tqs is even lower and the duration of the QS regime 
is longer, in some way mimicking larger systems |13|. We 
also performed simulations for equilibrium (EQ) initial 
conditions: waiting an appropriate transient after picking 
spatial coordinates and momenta from Boltzmann-Gibbs 
statistics. 

Integration of Eqs. J5J) yields phases in (—00,00). 
Then, we measured the histograms of phases at differ- 
ent times, accumulating the data over several realiza- 
tions to improve the statistics. Actually, the dynam- 
ics depends on the phases modulo 2tt, since they enter 
into the equations of motion as argument of sine and 
cosine functions. However, the statistics of unbounded 
phases is relevant as it reflects features of momentum 
space, which displays anomalies such as non-Maxwellian 
distributions 5] and two-time correlations signaling "ag- 
ing" effects For instance, a simple mathematical rela- 
tionship exists between the standard deviation of angles 
and velocity correlations [To). 

For low energies, phases are confined, while for suffi- 
ciently large energies, they evolve in diffusive motion. In 
Fig. n we show the properly scaled probability density 
functions (PDFs) of rotor phases at different times, for 
fixed size N — 1000 and energy per particle e = 0.69, 
starting from water-bag initial conditions. After tran- 
sient stages, numerical histograms can be very well de- 
scribed in the whole range by the g-Gaussian function Q 

P(6) = A{l + (q-l)(e/f3) 2 )^ , (3) 

where A is a normalization factor and (3 a positive con- 
stant. This function includes the standard Gaussian 
when q — > 1 and presents power-law tails for q > 1. Re- 
calling that the PDF © has variance a 2 = /? 2 /(5 - 3g), 
for q < 5/3, and considering normalized angles <p = 
then Eq. J3J can be written as the uni-parametric func- 
tion 

P^)=A g (l + !fl/)"\ (4) 

with A i = y/ *(5-3q) r(i7y-i7-i/2) ■ At each instant 1 01 
the dynamics, a 2 is computed as 

o 2 {t) = {(0- {0)t) 2 )u (5) 

where (...}* denotes average over the N rotors and over 
different realizations of the dynamics at time t. As far 
as we know, this is the first time that pure g-Gaussian 
PDFs arising from a Hamiltonian dynamics are detected. 

In Fig. [2Ji, the evolution of q for different N is dis- 
played, including the fitting values used in Fig. ^ for 
N = 1000. Parameter q increases up to a steady value in 



the long-time limit, that for all N falls within the range 
q ~ 1.51 ± 0.02. In the same figure, we also exhibit 
the temporal evolution of temperature T(t) — J2i(Pi)/N 
(Fig. |2| a) and dispersion a(t) (Fig. |2Jb), in order to es- 
tablish a parallel between the different stages of the rel- 
evant quantities. Instead of single runs, averaged quan- 
tities are presented. Since they are much less noisy, they 
allow the representation of several curves in the same plot 
without distorting the mean features observed in individ- 
ual runs. We can see that q attains a steady value ap- 
proximately when the QS— >EQ transition is completed. 
In Fig.|31the same data are presented for rescaled times. 

The PDFs of rotor phases, for N = 500, and e = 0.69, 
but starting from an equilibrium configuration are shown 
in Fig. 01 In this case, histograms present pronounced 
shoulders that persist for long times and can not be well 
described by g-Gaussian functions. However, as times 
goes by, the shoulders shift away from the center and the 
histograms tend to a g-Gaussian function, with g going 
to q ~ 3/2 from above in the long-time limit. 

Then, starting both from equilibrium and out-of equi- 
librium (water-bag) configurations, initially confined 
densities of phases develop power-law tails and adopt q- 
Gaussian forms. 

Diffusion of spatial coordinates may be characterized 
by the average squared displacement a 2 (t) of the angles 9 
as defined by Eq. JSJ. In the one-dimensional generalized 
Einstein relation 

a 2 (i) = ID t** , (6) 
where D is the diffusion constant, the case 7 = 1 corre- 
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FIG. 1: Histograms of rotor phases at different instants of 
the dynamics (symbols). Simulations for N = 1000 were per- 
formed starting from regular water-bag initial conditions at 
e = 0.69 (conditions leading to QS states). Countings were 
accumulated over 100 realizations, at times tk = 2 fc , with 
k — 6, 8, ..14, growing in the direction of the arrow up to 
t = 16384. Solid lines correspond to g-Gaussian fittings. His- 
tograms were shifted for visualization. Inset: log-log repre- 
sentation of the fitted data. 
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FIG. 3: Averaged time series of (a) temperature T, (b) lo- 
cal exponent 7 and parameter q (symbols), as a function of 
t/N 1 - 7 . Data are the same presented in Fig. |2] 



FIG. 2: Averaged time series of (a) temperature T, (b) de- 
viation a, (c) diffusion exponent 7 and (d) parameter q, for 
e = 0.69 and different values of N (N = 500 x 2 fe , with 
fc = 0, ...,9). Bold lines correspond to N — 500, as refer- 
ence, and iV increases in the direction of the arrows up to 
N = 256000. Averages were taken over 2.56 x 10 5 /N real- 
izations, starting from a waterbag configuration at t — 0. In 
panel (d), the fitting error is approx. 0.03. Dotted lines are 
drawn as references. In (a), they correspond to temperatures 
at equilibrium (Teq = 0.476) and at QS states in the TL 
(Tqs = 0.38). In (b), to ballistic motion (7 = 2) and normal 
diffusion (7 = I). 




sponds to normal diffusion, 7 < 1 to sub-diffusion and 
super-diffusion occurs for 7 > 1. The evolution of a is 
shown in Figs. [3) and|Sj5, for water-bag and equilibrium 
initial preparations, respectively. In order to detect dif- 
ferent regimes, it is useful to obtain an instantaneous 
exponent 7 as a function of time by taking the logarithm 
in both sides of Eq. (Jfjjl and differentiating with respect 
to hit: 



FIG. 4: Histograms of rotor phases at different instants of the 
dynamics (symbols). Simulations were performed for N = 500 
and e = 0.69, starting from an equilibrated initial condition. 
Countings were accumulated over 200 realizations, at times 
t k = 0.1 x 4 fc , k = 3,4, ..10, growing in the direction of the 
arrow up to t ~ 1.05 x I0 5 ). The g-Gaussian function with 
q = 1.53 was plotted for comparison (solid line). Histograms 
were shifted for visualization. 



7 (t) = 



d(lna 2 ) 
rf(lnt) 



(7) 



The outcome of this procedure can be seen in Fig. 
for water bag initial conditions. The same analysis for 
systems prepared in an equilibrated configuration is pre- 
sented in Fig. [SJ:. 

For the latter class of initial conditions, phase motion 



is ballistic (7 = 2) at short time scales in which rotors 
move almost freely, while phases display normal diffusion 
(7 = 1) in the long-time limit, as shown in Fig. [5] (see 
also [TJ| ). The crossover between both behaviors shifts 
to larger times as N increases. We observed this same 
pattern of behaviors also at supercritical energies (e = 5), 
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FIG. 5: Time series of (a) temperature T, (b) deviation a and 
(c) diffusion exponent 7, for N = 500 and e = 0.69. At t = 0, 
the system is in an equilibrium configuration. The outcome 
of a single run is exhibited. 



although in Ref. 0] only the ballistic regime was reported 
since normal diffusion settles at times longer than those 
analyzed in that work. 

For water-bag initial conditions, both ballistic and nor- 
mal regimes are also observed for short and long times, 
respectively. However, in this case, the crossover is more 
complex (see Figs. Hfc) and Etc)): Exponent 7 varies 
non-monotonously taking super-diffusive values and it 
does not present a well defined plateau as temperature 
does. For this intermediate stage, super-diffusion has 
been reported before and attributed to a kind of Levy- 
walk mechanism, yielding a succession of free walks and 
trapping events Also, it has been interpreted from a 
topological perspective (lq. 

The main features one observes in the profile of local 
exponent 7 vs. time can be summarized as follows: (i) In 
a first stage, 7 takes a maximum value that remains close 
to 7 = 2, corresponding to ballistic motion. This stage 
occurs at the beginning of the QS regime and its duration 
increases with N. (ii) Then 7 reaches a minimum value 
that as N increases falls within the QS range. Also as N 
increases, the valley broadens and flattens, and the min- 
imum value tends to unity (see Fig.[3J). Thus, anomalous 
diffusion in QS states is a finite-size effect, as already 
pointed out by Yamaguchi from a different perspective 
and for a different initial preparation!^, (iii) Another 
maximum appears in correspondence to the rapid relax- 



ation from QS to EQ. The maximal value grows with N, 
largely exceeding the value 7 = 2. But the peak narrows, 
hence it lasts less. This peak corresponds approximately 
to the inflexion point in the time evolution of temperature 
(Figs. j2Jand[SJl, whose slope increases with N. In fact, in 
the QS^EQ relaxation, the rapid increase of T (kinetic 
energy) leads to an accelerated increase of the phases in 
average, (iv) In the final stage, 7 relaxes asymptotically 
to unity, indicating normal diffusion at very long times. 

III. DISCUSSION AND REMARKS 

From the analysis of 7, one concludes that diffusive mo- 
tion in QS trajectories strongly depends on the system 
size. In the 7 vs. t profile, as N increases, the neigh- 
borhood of the first minimum (that becomes contained 
in the QS range defined by T) flattens, defining a quasi- 
stationary value of 7. This steady value goes to one in 
the TL. All these features are neatly observed in Fig. [3] 
Then, diffusion in QS states becomes of the normal type 
in that limit. For the whole time span, non-trivial val- 
ues of 7 occur in the TL only at non-stationary (nor 
quasi-stationary) stages, consistently with observations 
by Yamaguchi |lOj. 

Moreover, the spreading of initially confined phases de- 
velops power-law tails. For QS trajectories, histograms 
adopt a q-Gaussian form with q increasing up to a steady 
value q ~ 3/2. A similar steady value is also attained 
asymptotically when starting at equilibrium. For both 
classes of initial conditions, stabilization of parameter q 
is reached only when equilibrium holds and 7 attains the 
steady value corresponding to normal diffusion. 

Let us note that Levy densities present power-law tails 
with exponents restricted to a given interval (yielding 
divergent second moment) that in terms of parameter 
q corresponds to q > 5/3, therefore, not including the 
observed values. While the generalization of the stan- 
dard diffusion equation with fractional spatial derivatives 
leads to Levy functions, the non-linear generalization 
d t P{x,t) = Dd xx [P(x, t)] 2 ~ q , with constant D and q, 
has <?-Gaussian functions as long-time solutions pl|. The 
non-linear equation yields anomalous diffusion of the cor- 
related type with exponent 7 = 2/(3 — q). In our case, 
power-law tails also develop. Although they are charac- 
terized by q changing with time, the fact that it attains 
a steady value suggests that the spreading of phases may 
be ruled by an akin non-linear process. In such case, 
effective parameters are q e ff ~ 1.5 and 7 e // — 1.33. 
However, this issue as well as a possible connection with 
Nonextensive Statistics should be further investigated. 
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